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Description 

BIOLOGICAL DATA SET COMPARISON METHOD 
Technical Field 

The technical field relates to methods of identifying common properties 
within a set of biomolecules and properties that connect two or more sets of 
biomolecules, and also relates to methods for deriving functional explanations 
or hypotheses to explain the relationship between a set of biomolecules (e.g., 
genes, proteins) and between multiple sets of biomolecules. 

Table of Abbreviations 



3D three-dimensional 

BIOS basic input/output system 

BLAST Basic Local Alignment Search Tool 

CGI common gateway interface 

cM centimorgan 

DNA deoxyribonucleic acid 

HSPs high scoring sequence pairs 

LAN local area network 

LOD Log of the odds ratio 

NCBI National Center for Biotechnology 
Information 

NLM National Library of Medicine 

PGR polymerase chain reaction 

PNA peptide nucleic acid 

OMIM Online Mendelian Inheritance in Man 

RAM random access memory 

rmsd root-mean-squared distance 

RNA ribonucleic acid 

ROM read only memory 

SAN system area network 

URL uniform resource locator 
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USB universal serial bus 

WAN wide area network 



Amino Acid Abbreviations and Corresponding mRNA Codons 



Amino Acid 


3-Letter 


1 -Letter 


mRNA Codons 


Alanine 


Ala 


A 


GCA GCC GCG GCU 


Arginine 


Arg 


R 


AGA AGG CGA CGC CGG CGU 


Asparagine 


Asn 


N 


AAC AAU 


Aspartic Acid 


Asp 


D 


GAC GAU 


Cysteine 


Cys 


C 


UGC UGU 


Glutamic Acid 


Glu 


E 


GAA GAG 


Glutamine 


Gin 


Q 


CAACAG 


Glycine 


Gly 


G 


GGA GGC GGG GGU 


Histidine 


His 


H 


CAC CAU 


Isoleucine 


He 


I 


A 1 1 A All A 1 III 

AUA AUC AUU 


Leucine 


Leu 


L 


UUA UUG CUA CUC CUG CUU 


Lysine 


Lys 


K 


AAA AAG 


Methionine 


Met 


M 


AUG 


Proline 


Pro 


P 


CCA CCC CCG CCU 


Phenylalanine 


Phe 


F 


UUC UUU 


Serine 


Ser 


S 


ACG AGU UCA UCC UCG UCU 


Threonine 


Thr 


T 


ACAACC ACGACU 


Tryptophan 


Trp 


W 


UGG 


Tyrosine 


Tyr 


Y 


UAC UAU 


Valine 


Val 


V 


GUA GUC GUG GUU 



Background Art 

Biomedical research is in the midst of an unprecedented data 
explosion. Complete genome sequences of prokaryotic organisms are 
30 appearing in the literature and on the World Wide Web on almost a weekly 
basis. See e.g., http://igweb.integratedgenomics.com/GOLD/. Several 
complete genomes from model eukaryotic organisms have also been 
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sequenced, and many more sequencing projects are in various stages of 
planning and execution See e.g., http://www.nih.gov/science/models/. The 
sequence of the human genome is also now freely available in "finished" form. 
See e.g., http://www.ncbi.nlm.nih.gov/genome/guide/human/ or 
5 http://www.ensembl.org/Homo_sapiens/. Combined with the growing 
availability of high-throughput and genome-wide experimental methods, this 
deluge of data facilitates the potential for comparisons of sequence, structure, 
mRNA- or protein-expression levels, and function between all human genes 
and the genes of model organisms. It also opens up new challenges for 

10 determining the functional and cellular role for the many as yet 
uncharacterized genes within these organisms. 

As research into genomics and proteomics progresses, experimental 
results are beginning to transcend a single gene of interest and are more 
commonly involving sets of genes or other biomolecules that behave in some 

15 sense "similarly" or share a common property. Although computational tools 
that allow for a comparison of one gene to all other known genes at the level 
of primary nucleic acid or amino acid sequence have existed for some time 
(e.g., BLAST; Altschul ef a/., 1990), such comparisons often do not yield 
sufficient information to allow for the identification of a specific function for that 

20 gene. Indeed, it is very common for genes that share little or no similarity at 
the nucleic acid sequence level to encode proteins that have related functions 
or roles. For example, two genes might encode enzymes that catalyze 
adjacent steps in the same biochemical pathway, and the functional disruption 
of either gene might lead to a similar outcome for the cell or organism (e.g., a 

25 human disease). These genes would be unlikely to exhibit similarity at the 
primary nucleic acid sequence level, and thus current search strategies would 
not identify these genes as being related despite the similar phenotype that 
would result from their functional disruption. By way of additional example, 
this problem is also encountered in areas such as transcriptome analysis, 

30 where lists of genes with similar expression levels or time-profiles are 
generated from each experiment. Thus, there persists a great need for 
computational methods for determining the underlying commonality among a 
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set of genes and for ways of assigning consensus annotations to such gene 
sets. 

One currently available approach for analyzing genes is a World Wide 
Web-based tool that collects and displays information gene-by-gene for a 
5 predefined set of genes, such as disease candidates by creating a "home 
page" for each gene in the set. Halushka et a/., 1999. This and other 
approaches (see e.g., Bouton and Pevsner, 2000; Bouton and Pevsner, 2002; 
Khatri et a/., 2002; Ostermeier et a/, 2002) lack breadth and do not 
comprehensively address the universe of possible interactions, traits, and 

1 0 characteristics between genes. 

Some other approaches involving text mining of published scientific 
abstracts have been developed for use in gene expression profiling (see e.g., 
Tanabe ef a/., 1999; Masys et a/., 2001; Blaschke ef a/ M 2001), or for finding 
links between genes and diseases (Jenssen et a/., 2001; Perez-lratxeta ef a/., 

15 2002a). The latter group has recently demonstrated the feasibility of mining 
MEDLINE abstracts to generate lists of candidate genes that are believed to 
be associated with a group of inherited diseases. Perez-lratxeta et a/., 2002b. 

Computational methods have been proposed that pertain to partitioning 
of genotype variation into clusters that predict quantitative trait variation, such 

20 as elevated plasma triglyceride levels. Nelson et a/., 2001. An extension of 
this method has been used to uncover a combination of polymorphisms in 
several estrogen metabolism genes that correlates with increased sporadic 
breast cancer occurrence. Ritchie et a/., 2001. A support-vector machine 
approach was employed to make gene functional classifications based on 

25 phylogenetic profiles and expression data. Pavlidis et a/., 2002. Additionally, 
a graph theoretic method for combining microarray and data with protein 
interaction maps as a way of annotating sets of genes from transcriptome 
experiments has been described, del Rio et a/., 2001 . 

While the above methods attempt to address the general problem of 

30 assigning consensus annotations to gene sets, these approaches do not offer 
a comprehensive solution to the problem of identifying the properties of a set 
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of biomolecules and correlating these properties with other sets of 
biomolecules for which a common property has been defined. 

What is needed, therefore, is a method of identifying various properties 
of a given set of biomolecules and correlating these properties with multiple 
5 sets of biomolecules that are common to a given biological process or 
pathway. Such a method would facilitate the characterization of a set of 
unknown biomolecules, including an assessment of the function of the 
unknown biomolecules. These and other problems are addressed herein. 

10 Summary 

Provided is a method of identifying a relationship between one or more 
candidate biomolecules and one or more reference biomolecules. In one 
embodiment, the method comprises: (a) inputting to a computer a query set 
describing the one or more candidate biomolecules; (b) comparing the query 

15 set with a target database describing the one or more reference biomolecules, 
wherein the one or more reference biomolecules are grouped into one or 
more buckets, and wherein the one or more reference biomolecules of each 
bucket share a common property; (c) counting a number of matches between 
each query set and each bucket of the target database; and (d) statistically 

20 analyzing each match, wherein the presence of a statistically significant match 
identifies a relationship between the query set and a bucket of the target 
database. 

Also provided is a method of identifying a relationship between two or 
more region sets, each region set describing one or more candidate 

25 biomolecules, and a target database describing one or more reference 
biomolecules grouped into one or more buckets. In one embodiment, the 
method comprises: (a) providing a query set describing two or more 
region sets, each region set comprising one or more candidate biomolecule 
sequences extracted from one genetic region; (b) comparing the query set 

30 with target database sequences describing one or more reference 
biomolecule sequences, wherein the target database sequences grouped into 
one or more buckets, and wherein the one or more reference biomolecules of 
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each bucket share a common property; (c) counting a number of matches 
between each query set and each bucket of the target database; and (d) 
statistically analyzing each match, wherein the presence of a statistically 
significant match identifies a relationship between the query set and a bucket 
5 of the target database. In one embodiment, the method further comprises (e) 
constructing a plurality of replicates of the one or more query sets; (f) 
modeling the replicates at random chromosomal locations to form a random 
location data set; (g) processing the random location data set by following 
steps (a) - (d); (h) quantifying the number of times each match is found to 
10 surpass a predetermined threshold to form a statistically significant set of 
random location matches; and (i) comparing the statistically significant set of 
random location matches to the statistically significant relationship of steps (a) 
-(d). 

In various embodiments, query sets comprise one or more sequences, 

15 including, but not limited to, DNA, RNA, or protein sequences. In one 
embodiment, these sequences are derived from one genetic region. In one 
embodiment, the one or more candidate biomolecules and the one or more 
reference biomolecules are all selected from the group consisting of proteins, 
nucleic acids, and small molecules. In one embodiment, the comparing 

20 comprises employing a BLAST-based algorithm to identify similar or identical 
sequences. In one embodiment, the counting comprises applying one or 
more principles chosen from the group consisting of (a) each query set 
candidate sequence can match at most one reference sequence in any given 
bucket; (b) each query set candidate sequence can possess a match in one 

25 or more different buckets; and (c) once a candidate sequence in the query set 
matches a specific bucket reference sequence in the target database, any 
subsequent matches of that same candidate sequence to other reference 
sequences in that bucket do not increase the match count for the bucket. In 
one embodiment, the statistically analyzing comprises computing one or more 

30 statistics for each match, which can optionally be sorted and/or outputted to a 
webpage comprising one or more hyperlinks. 
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Also provided is a computer-readable medium having stored thereon a 
data structure having multiple data fields, comprising (a) a first data field 
containing data representing a bucket; (b) a second data field containing data 
representing a name for the bucket; and (c) a third data field containing data 
5 representing a list of members of the bucket, wherein the members have a 
common property. 

Also provided is a method of making a target database. In one 
embodiment, the method comprises: (a) identifying a source of informative 
content; (b) arranging informative content from the source of informative 

10 content into a set of buckets, wherein the buckets are given names; (c) 
gathering the names of the buckets and a list of biomolecules present in each 
bucket; and (d) creating and loading into a database data fields containing 
data representing (i) the set of buckets; (ii) the list of biomolecules present in 
each bucket; and (iii) a description for each biomolecule present in each 

15 bucket. In one embodiment, the source of informative content is a publicly 
available database, including, but not limited to, SwissProt, TrEMBL, and 
NCBI. In one embodiment, the gathering is accomplished using a source- 
specific parsing script. In one embodiment, the creating and loading is 
accomplished using a database loading script. In one embodiment, the data 

20 representing a description for each biomolecule present in each bucket is 
selected from the group consisting of a nucleic acid sequence, an amino acid 
sequence, or an identification number, wherein the identification number 
allows for retrieval of a nucleic acid sequence or an amino acid sequence. 

Also provided is a computer readable storage device embodying 

25 programs of instructions executable by a computer for performing the 
disclosed methods. 

Accordingly, it is an object to provide a novel method for characterizing 
a set of biomolecules. This and other objects are achieved in whole or in part 
as disclosed herein. 

30 An object having been stated hereinabove, other objects will be evident 

as the description proceeds, when taken in connection with the accompanying 
drawings and examples as best described hereinbelow. 
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Brief Description of the Drawings 
Figure 1 illustrates an exemplary general purpose computing platform 
100 upon which the methods and systems disclosed herein can be 
5 implemented. 

Figure 2 is a flowchart of a process 200 for implementing the methods 
disclosed herein. 

Figure 3 is a flowchart of a process 300 for implementing a method of 
identifying a relationship between two or more regions sets. 
10 Figure 4 is a database relation diagram 400 showing exemplary data 

that is stored in each field and how the data in one field relates to the data in 
another field. 

Detailed Description 

15 The disclosed methods and data structures can be implemented in 

hardware, firmware, software, or any combination thereof. In one exemplary 
embodiment, the methods and data structures disclosed herein for classifying 
biomolecules can be implemented as computer readable instructions and data 
structures embodied in a computer-readable medium. 

20 With reference to Figure 1 , an exemplary system includes a general 

purpose computing device in the form of a conventional personal computer 
100, including a processing unit 101, a system memory 102, and a system 
bus 103 that couples various system components including the system 
memory to the processing unit 101. System bus 103 can be any of several 

25 types of bus structures including a memory bus or memory controller, a 
peripheral bus, and a local bus using any of a variety of bus architectures. 
The system memory includes read only memory (ROM) 104 and random 
access memory (RAM) 105. A basic input/output system (BIOS) 106, 
containing the basic routines that help to transfer information between 

30 elements within personal computer 100, such as during start-up, is stored in 
ROM 104. Personal computer 100 further includes a hard disk drive 107 for 
reading from and writing to a hard disk (not shown), a magnetic disk drive 108 
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for reading from or writing to a removable magnetic disk 109, and an optical 
disk drive 110 for reading from or writing to a removable optical disk 111 such 
as a CD ROM or other optical media. 

Hard disk drive 107, magnetic disk drive 108, and optical disk drive 1 10 
5 are connected to system bus 103 by a hard disk drive interface 112, a 
magnetic disk drive interface 113, and an optical disk drive interface 114, 
respectively. The drives and their associated computer-readable media 
provide nonvolatile storage of computer readable instructions, data structures, 
program modules, and other data for personal computer 100. Although the 

10 exemplary environment described herein employs a hard disk, a removable 
magnetic disk 109, and a removable optical disk 111, it will be appreciated by 
those skilled in the art that other types of computer readable media which can 
store data that is accessible by a computer, such as magnetic cassettes, flash 
memory cards, digital video disks, Bernoulli cartridges, random access 

15 memories, read only memories, and the like can also be used in the 
exemplary operating environment 

A number of program modules can be stored on the hard disk, 
magnetic disk 109, optical disk 111, ROM 104, or RAM 105, including an 
operating system 115, one or more applications programs 116, other program 

20 modules 117, and program data 118. 

A user can enter commands and information into personal computer 
100 through input devices such as a keyboard 120 and a pointing device 122. 
Other input devices (not shown) can include a microphone, touch panel, 
joystick, game pad, satellite dish, scanner, or the like. These and other input 

25 devices are often connected to processing unit 101 through a serial port 
interface 126 that is coupled to the system bus, but can be connected by 
other interfaces, such as a parallel port, game port or a universal serial bus 
(USB). A monitor 127 or other type of display device is also connected to 
system bus 103 via an interface, such as a video adapter 128. In addition to 

30 the monitor, personal computers typically include other peripheral output 
devices, not shown, such as speakers and printers. The user can use one of 
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the input devices to input data indicating the user's preference between 
alternatives presented to the user via monitor 127. 

Personal computer 100 can operate in a networked environment using 
logical connections to one or more remote computers, such as a remote 
5 computer 129. Remote computer 129 can be another personal computer, a 
server, a router, a network PC, a peer device or other common network node, 
and typically includes many or all of the elements described above relative to 
personal computer 100, although only a memory storage device 130 has been 
illustrated in Fig. 1. The logical connections depicted in Fig. 1 include a local 
10 area network (LAN) 131, a wide area network (WAN) 132, and a system area 
network (SAN) 133. Local- and wide-area networking environments are 
commonplace in offices, enterprise-wide computer networks, intranets and the 
Internet. 

System area networking environments are used to interconnect nodes 

15 within a distributed computing system, such as a cluster. For example, in the 
illustrated embodiment, personal computer 100 can comprise a first node in a 
cluster and remote computer 129 can comprise a second node in the cluster. 
In such an environment, it is preferable that personal computer 100 and 
remote computer 129 be under a common administrative domain. Thus, 

20 although computer 129 is labeled "remote", computer 129 can be in close 
physical proximity to personal computer 100. 

When used in a LAN or SAN networking environment, personal 
computer 100 is connected to local network 131 or system network 133 
through network interface adapters 134 and 134a. Network interface 

25 adapters 134 and 134a can include processing units 135 and 135a and one or 
more memory units 136 and 136a. 

When used in a WAN networking environment, personal computer 100 
typically includes a modem 138 or other device for establishing 
communications over WAN 132. Modem 138, which can be internal or 

30 external, is connected to system bus 103 via serial port interface 126. In a 
networked environment, program modules depicted relative to personal 
computer 100, or portions thereof, can 1 be stored in the remote memory 
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storage device. It will be appreciated that the network connections shown are 
exemplary and other approaches to establishing a communications link 
between the computers can be used. 

5 L Definitions 

Following long-standing patent law convention, the terms "a M and "an" 
mean "one or more" when used in this application, including the claims. 

As used herein, the term "about," when referring to a value or to an 
amount of mass, weight, time, volume, concentration or percentage is meant 

10 to encompass variations of ±20% or ±10%, in another example ±5%, in 
another example ±1%, and in still another example ±0.1% from the specified 
amount, as such variations are appropriate to perform the disclosed method. 

As used herein, the terms "amino acid" and "amino acid residue" are 
used interchangeably and mean any of the twenty naturally occurring amino 

15 acids. An amino acid is formed upon chemical digestion (hydrolysis) of a 
polypeptide at its peptide linkages. In keeping with standard polypeptide 
nomenclature, abbreviations for amino acid residues are shown in tabular 
form presented hereinabove. In addition, the phrases "amino acid" and 
"amino acid residue" are broadly defined to include modified and unusual 

20 amino acids. 

As used herein, the term "biomolecule" means any molecule isolated 
from, derived from, or based on a molecule found in a living organism, 
including viruses. The term biomolecule includes, but is not limited to, both 
proteins and nucleic acids (RNA and DNA). Biomolecules can be polymeric in 

25 nature and can comprise a unique sequence of monomers; for example, a 
biomolecule can comprise a nucleic acid (e.g., a gene, and fragments 
thereof), an amino acid, a derivatized protein (e.g., a glycosylated protein), a 
nucleic acid comprising a nucleic acid analog, a peptide nucleic acid (PNA), 
an antibody, as well as peptides, polypeptides, proteins and fragments 

30 thereof. As used herein, the term "biomolecule" also refers to any molecule 
that is capable of producing a biological effect or participating in a biological 



WO 2005/003308 



PCT/US2004/019932 



-12- 

process. In this context, a biomolecule includes, but is not limited to, a small 
molecule such as a drug. 

As used herein, the term "BLAST-formatted database" means a 
database wherein the data representing a nucleic acid or amino acid 

5 sequence of a candidate or reference biomolecule is in a form amenable to 
manipulation by BLAST and BLAST-based algorithms. The proper form for 
such sequences is described in Altschul et a/., (1990). See also 
http://blast.wustl.edu/doc/FAQ-lndexing.html. The BLAST-formatted database 
acts as a master repository for all nucleic acid and amino acid sequences. It 

10 includes data entries for nucleic acid and amino acid sequences 
corresponding to all reference biomolecules as well as identification or 
accession numbers by which these sequences can be accessed for use in the 
methods and devices disclosed herein. In addition, data is automatically 
added to the BLAST-formatted database corresponding to the nucleic acid 

15 and amino acid sequences of all candidate biomolecules. 

As used herein, the term "bucket" means any grouping of biomolecules 
(e.g., genes or gene products) that share a biological property. For example, 
a gene or gene product can have an identifier and/or an associated sequence 
(amino acid or nucleic acid). In one example, an identifier is a standard name 

20 for the gene or gene product (e.g., "human beta-globin"). In another example, 
the identifier is an identification number or an accession number that allows 
the sequence of the gene or gene product to be retrieved from a source (e.g., 
the NCBI accession number for the human beta-globin complete coding 
sequence is AF007546). A source includes, but is not limited to a public or 

25 private database. The identifier need not be unique, and a given gene can be 
a member of one or more buckets. 

Each bucket can have a unique name, which can also indicate its origin 
and/or creator. Buckets and collections of buckets can be created by 
individuals or they can be defined as the results from various types of 

30 analyses. For example, a bucket can comprise a set of genes found to be 
more highly expressed in a particular tumor cell compared with a normal cell. 
Buckets can also be created from public-domain databases. As an additional 
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example, a bucket can include all the component enzymes in a metabolic 
pathway, all the protein components in a signaling pathway, biomolecules 
mentioned in the same publication, biomolecules mentioned in publications on 
the same subject, sets of proteins sharing a particular sequence motif or 
5 domain, sets of genes known to be present on an oligonucleotide array or 
chip, genes classified into particular categories according to an ontology, 
gene products present in a particular tissue or organ or subcellular location, or 
genes in which a particular keyword occurs somewhere in their associated 
annotations. A bucket can form an element of a target database. 

10 As used herein, the term "bucket source" means any medium or entity 

to which the origin of the bucket can be traced. For example, a bucket source 
can be a user. In another example, a bucket source can be a database. In 
yet another example, a bucket source can be the results of a search of a 
database done with user-specified parameters. Defining a bucket source can 

15 be useful as an approach for identifying different buckets that have the same 
name. The use of bucket sources also allows broad categories of buckets to 
be defined, such as "pathway" or "function" buckets. 

As used herein, the terms "candidate biomolecule" and "candidate 
sequence" are used interchangeably, and mean a biomolecule or sequence 

20 that is part of a query set to be compared to a target database. Candidate 
biomolecules are ones that the user is attempting to characterize as having or 
not having the various properties that are represented by the buckets of the 
target database. This characterization is accomplished by comparing a 
candidate biomolecule to the reference biomolecules of the target database 

25 and statistically analyzing the number of matches that result from the 
comparison. When a statistically significant match (of the query set) is found 
to a particular bucket, the user can infer that the candidate biomolecule has 
the property that is common to the reference biomolecules that are members 
of the bucket to which the match was made. 

30 As used herein, the terms "describing" and "description" as they relate 

to biomolecules mean any categorization of the biomolecule that relates to its 
identity or to a property it possesses. In one example, a biomolecule can be 
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described by its common name, such as "human beta-globin", "mouse 
erythropoietin receptor, "Drosophila fushi tarazu", etc. In another example, a 
biomolecule can be described by its nucleic acid or amino acid sequence. In 
another example, a biomolecule can be described by an identification number 

5 or accession number that allows its corresponding nucleic acid and/or amino 
acid sequence to be retrieved from a source such as a public or private 
database. In yet another example, a biomolecule can be described by a 
property that it possesses. In one example, a property can be a functional 
description of the biomolecule such as "kinase", "receptor", "cytokine", 

10 "oncogene", "ligand", etc. In another example, a property can include the 
organism from which the biomolecule was isolated. In another example, the 
property can include a biochemical pathway in which the gene product plays a 
role including, but not limited to pyrimidine biosynthesis, the citric acid cycle, 
fatty acid biosynthesis, the pentose cycle, amino acid biosynthesis, etc. In yet 

15 another example, the property can include a three-dimensional (3D) structural 
feature of the biomolecule. Several ways of incorporating structural 
information into a search exist including, but not limited to creating sets of 
buckets based on public databases that impose a structural hierarchy on 
those proteins with known three-dimensional structures. Exemplary public 

20 databases include CATH (http://www.biochem.ucl.ac.uk/bsm/cath_new/index. 
html) and SCOP (http://scop.berkeley.edu/). It is generally believed that 
proteins with at least 30% overall amino acid identity are likely to fold into very 
similar structural conformations (McGuffin & Jones 2002). 

A more general method might be to reduce or project known 3D 

25 structures to a sequence-like character string, comprising the secondary 
structure adopted by each amino acid (e.g., hhhhhhhhhsssshhhhhhhhhhhhh 
as a helix-loop-helix motif). A BLAST-like method could optionally be used to 
compare the length and order of secondary-structural elements of known 
proteins. Secondary structure predictions for proteins with no known structure 

30 could also be compared to those of a database of known structures (see 
Aurora & Rose 1998). 

Yet another possibility is to create structure-specific buckets by 
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computing a root-mean-squared distance (rmsd) measure between the 3D 
structural coordinates of any two proteins. For example, buckets for all 
structures within 2 A rmsd of each other could be defined. 

As used herein, the phrase "extracted from one genetic region" refers 

5 to sequences derived from genes that are present in a contiguous region of a 
genome or to protein sequences that are encoded by sequences derived from 
genes that are present in a contiguous region of a genome. "One genetic 
region" and "the same region of a genome" include, but are not limited to a 
chromosome, an arm of a chromosome, a portion of a chromosome contained 

10 between two markers, and a band of a chromosome as visualized by banding 
techniques that are known in the art such as Giemsa banding. These terms 
also include any other measure of physical proximity on a chromosome, 
including but not limited to a kilobase, a megabase, or a centimorgan (cM). 

As used herein, the term "mutation" carries its traditional connotation 

15 and means a change, inherited, naturally occurring, or introduced, in a nucleic 
acid or polypeptide sequence, and is used in its sense as generally known to 
those of skill in the art. A mutation can be any (or a combination of) 
detectable, unnatural change affecting the chemical or physical constitution, 
mutability, replication, phenotypic function, or recombination of one or more 

20 deoxyribonucleotides. Nucleotides can be added, deleted, substituted for, 
inverted, or transposed to new positions with and without inversion. Mutations 
can occur spontaneously and can be induced experimentally by application of 
mutagens. A mutant variation of a nucleic acid molecule results from a 
mutation. A mutant polypeptide can result from a mutant nucleic acid 

25 molecule and can also refer to a polypeptide that is modified at one or more 
amino acid residues from the wild-type (/.e., naturally occurring) polypeptide. 
For example, the mutation can be a point mutation or the addition, deletion, 
insertion, and/or substitution of one or more nucleotides, or any combination 
thereof. The mutation can be a missense or frameshift mutation. 

30 Modifications can be, for example, conserved or non-conserved, natural or 
unnatural. 

As used herein, "nucleic acid" and "nucleic acid molecule" refer to any 



WO 2005/003308 



PCT/US2004/0 19932 



-16- 

of deoxyribonucleic acid (DNA), ribonucleic acid (RNA), oligonucleotides, 
fragments generated by the polymerase chain reaction (PCR), and fragments 
generated by any of ligation, scission, endonuclease action, and exonuclease 
action. Nucleic acids can comprise monomers that are naturally occurring 
5 nucleotides (such as deoxyribonucleotides and ribonucleotides), or analogs of 
naturally occurring nucleotides (e.g., a-enantiomeric forms of naturally 
occurring nucleotides), or a combination of both. Modified nucleotides can 
have modifications in sugar moieties and/or in pyrimidine or purine base 
moieties. Sugar modifications include, for example, replacement of one or 

10 more hydroxyl groups with halogens, alkyl groups, amines, and azido groups. 
Sugars can also be functionalized as ethers or esters. Moreover, the entire 
sugar moiety can be replaced with sterically and electronically similar 
structures, such as aza-sugars and carbocyclic sugar analogs. Examples of 
modifications in a base moiety include alkylated purines and pyrimidines, 

15 acylated purines or pyrimidines, or other well-known heterocyclic substitutes. 
Nucleic acid monomers can be linked by phosphodiester bonds or analogs of 
phosphodiester bonds. Analogs of phosphodiester linkages include 
phosphorothioate, phosphorodithioate, phosphoroselenoate, 

phosphorodiselenoate, phosphoroanilothioate, phosphoranilidate, 

20 phosphoramidate, and the like. The term "nucleic acid" also includes so- 
called "peptide nucleic acids," which comprise naturally occurring or modified 
nucleic acid bases attached to a polyamide backbone. Nucleic acids can be 
either single stranded or double stranded. 

As used herein, the term "property" denotes any feature of a 

25 biomolecule. Properties include, but are not limited to, sequence similarity 
and/or identity, chromosomal location, involvement in a particular biochemical 
pathway, association with genetic disease, expression in a context, three- 
dimensional structural features, and having or encoding a particular functional 
domain. Representative functional domains include, but are not limited to, 

30 kinase domains, growth factor binding domains, phosphorylation sites, 
glycosylation sites, protein and/or nucleic acid binding sites, protein-protein 
interaction domains, and post-translational modification sites. 
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As used herein, the term "quality checking" means the application of 
subjective criteria to assess the usefulness of a bucket. Quality checking 
ensures that all reference biomolecules that have been grouped into a bucket 
share the common property used to describe the bucket. These criteria 

5 attempt to take into account the nature of the data analysis involved in 
assembling the bucket. For example, reliable human-annotated sources (e.g., 
the SwissProt database) would receive a higher rank than a set generated by 
some automated computational procedure. 

As used herein, the term "query set" means any item or group of items 

10 arranged in such a way as to allow for comparison to a target database. By 
way of example and not limitation, a query set can include a nucleic acid 
sequence, an amino acid sequence, or a combination thereof. Query sets 
can be produced by manual grouping of items. In another example, a query 
set can be produced by techniques including, but not limited to text mining of 

15 sequence databases and literature, homology searches, annotation keyword 
searches, or any other technique that generates a group of items that are 
believed to share a common property. Query sets can comprise results from 
one or more biological experiments, for example as raw data or as a product 
of statistical or other data analyses. 

20 As used herein, the term "query sequence" means a member of a 

query set. In one example, a query sequence is a nucleic acid or amino acid 
sequence. In one embodiment, query sequences can be grouped together to 
form one or more query sets. The query set(s) is/are then compared to a 
target database that has been organized into buckets, the members of each 

25 bucket sharing a common property. 

As used herein, the term "reference biomolecule" refers to the 
members of the buckets that make up a target database. In one embodiment, 
a "reference biomolecule" is a "reference sequence". Reference biomolecules 
are arranged in a target database into buckets, wherein the reference 

30 biomolecules in each bucket share a common property. 

As used herein, the term "region set" or "regions sets" means a set 
containing at least some, and optionally all, of the known and predicted genes 
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that lie within a contiguous region of a genome. A region set might have as its 
members all the genes either known or predicted to reside in one example on 
a certain chromosome, in another example on one arm of a certain 
chromosome, in another example on a portion of a chromosome contained 

5 between two markers, in another example on that area of a certain 
chromosome corresponding to a particular chromosomal band as visualized 
by G-banding with Giemsa stain, or in yet another example within a certain 
number of basepairs of each other on a certain chromosome. The certain 
number of basepairs can be measured in bases, kilobases, megabases, or 

10 cM. 

As used herein, the term "relationship" means any association between 
one or more entities. Relationships include, but are not limited to nucleic acid 
and/or amino acid sequence similarity and/or identity, presence in the same 
region of a genome or being encoded by genes present in the same region of 
15 the genome, having the same or a similar function, containing or encoding a 
common functional domain, containing a common three-dimensional 
structural feature, association with a similar phenotype such as a disease 
state, involvement in the same biochemical pathway, and any combination 
thereof. 

20 As used herein, the term "relevant universe of all characterized 

sequences" means all sequences that have been characterized to an extent 
sufficient to allow the user to conclude that the corresponding biomolecules 
should or should not be placed into a bucket. This conclusion can be based 
upon an assessment or a hypothesis as to whether or not a given biomolecule 

25 has the property shared by the members of a given bucket. When several 
complementary or competing sources exist for assigning biomolecules to a 
bucket (e.g., kinase buckets as defined from several different sources or 
methods) all can be included rather than attempting to choose one "best" set. 
As used herein, the terms "significance" and "significant" relate to a 

30 statistical analysis of the probability that there is a non-random association, or 
a more unusual relationship, between two or more entities. In one example, 
"significance" refers to the probability that an observed relationship occurred 
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by chance. To determine whether or not a relationship is "significant" or has 
"significance", statistical manipulations of the data can be performed to 
calculate a probability, expressed as a "p-value". Those p-values that fall 
below a user-defined cutoff point are regarded as significant. In one example, 
5 a p-value less than or equal to 0.05, in another example less than 0.01 , in 
another example less than 0.005, and in yet another example less than 0.001 , 
are regarded as significant. 

The term "similarity" can be contrasted with the term "identity". 
Similarity is determined using an algorithm including, but not limited to, the 

10 BLAST-based algorithms or the GAP program (available from the University 
of Wisconsin Genetics Computer Group, now part of Accelrys Inc., San Diego, 
California, United States of America). "Identity", however, means a nucleic 
acid or amino acid sequence having the same nucleic acid or amino acid at 
the same relative position in a given family member of a gene family or in a 

15 homologous nucleic acid or amino acid in a different organism. Homology 
and similarity are generally viewed as broader terms than the term identity. 
Biochemically similar amino acids, for example, leucine/isoleucine or 
glutamate/aspartate, can be present at the same position in a biomolecule- 
these are not identical per se, but are biochemically "similar." These are 

20 referred to herein as conservative differences or conservative substitutions. 
This differs from a conservative substitution or mutation at the DNA level, 
which is defined as a change in a nucleic acid residue that does not result in a 
change in the amino acid codon encoded by the DNA at the altered position 
(e.g., TCC to TCA, both of which encode serine). 

25 As used herein, the term "size" as it relates to a query set, a target 

database bucket, a genome, or a relevant universe of all characterized 
sequences, means the number of members present in the referenced item. 
For example, the size of a query set or a target database bucket would be the 
number of candidate biomolecules or reference biomolecules that make up 

30 the query set or target database bucket, respectively. Similarly, the size of a 
genome is the number of genes present in a genome or the number of gene 
products encoded by those genes. Also similarly, the size of the relevant 
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universe of all characterized sequences is the number of sequences that have 
been characterized sufficiently such that a user can either include or exclude 
a given biomolecule from a given bucket based upon the biomolecule having 
or lacking the property shared by the members of the bucket. The "size of the 

5 relevant universe" will typically be less than or equal to the size of the 
genome. It is also possible to define an "effective size" for a bucket, or for an 
entire genome, by performing redundancy analysis. Thus, if several very 
closely related sequences exist within a bucket (several mutant versions of 
the same protein, for example), one can define the number of substantially 

10 different members to be the "effective size" for that bucket. A similar 
correction could be applied on a per-genome basis as well. 

As used herein, the term "source of informative content" means any 
source of information that describes a relationship between biomolecules or 
assigns a property to a biomolecule. A source of informative content includes, 

15 but is not limited to an annotated database of nucleic acid or amino acid 
sequences. In this example, the annotations can include references to 
suspected functions, expression patterns, homologs or orthologs from the 
same or different species, presence on a particular microarray chip or in a 
particular cDNA library, or presence on a particular chromosome or region of 

20 a chromosome. Other non-limiting sources of informative content include 
journal articles, public databases, web pages or trees, scientific abstracts 
and/or posters, technical data sheets, or personal communications. 
Experimental results, whether raw or resulting from prior analysis, can also be 
sources of informative content. 

25 As used herein, the term "target database" means a collection of 

descriptions of one or more reference biomolecules. The reference 
biomolecules described in the collection are arranged in the target database 
into one or more buckets, wherein the members of each bucket share a 
common property. The reference biomolecules are further arranged such that 

30 the members of a bucket can be compared to a query set. 
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IL Biomolecule Analysis 

A representative embodiment is adapted to identify properties that are 
common between a query set and a target database. The method can be 
employed, for example, to identify the function of a gene product of one or 
5 more genes that form a query set. 

Referring now to Figure 2, given a set of sequences {i.e., a query set) 
at steps ST202, ST204a, and ST204b comprising one or more query 
sequences, the query set is compared using the BLAST algorithms (Altschul 
et a/., 1990) to a target database comprising one or more sequences grouped 

10 into one or more buckets (Figure 2 at step ST206). Each member of a query 
set is compared to each member of a target database. In one example, an 
embodiment can be configured to define a stringent threshold for filtering 
sequence match results. As shown in step ST208 of Figure 2, in this 
configuration, if a query sequence possesses above-threshold matches to 

15 more than one sequence in the target database, only the best match in each 
bucket is counted and the count for that bucket is incremented by 1 as shown 
in Figure 2 at step ST210. However, due to redundancy of the target 
database, a matching sequence can be a member of several different buckets 
of the target database. The query set can also contain redundancies, so once 

20 a candidate query set sequence has been matched to a given reference 
target database sequence, any subsequent matches to that same reference 
target database sequence in that bucket are ignored as shown in Figure 2 at 
step ST212. Thus, irrespective of the size of the query set, a particular bucket 
can have no more matches than the number of reference sequences that the 

25 bucket contains. Once a candidate query set sequence is compared to all the 
reference target database sequences in a given bucket, that process is 
repeated for the candidate query set sequence with the reference target 
database sequences of the next bucket, as shown in Figure 2 at step ST214. 
Once a given candidate query set sequence has been compared to all 

30 reference target database sequences, the process is repeated for the next 
candidate query set sequence as shown in Figure 2 at step ST216. Once all 
candidate query set sequences have been compared to all reference target 
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database sequences, each bucket with a count greater than 1 is collected as 
shown in Figure 2 at step ST218. 

As shown in Figure 2 at step ST220, to account for the different sizes 
of the buckets (/.e. f the different numbers of members that each bucket 

5 contains), a hypergeometric-distribution statistic is computed to assess the 
significance of the results. In this manner, a query set that matches 49 of 50 
sequences in one bucket, for example, is considered to be a more significant 
result than a match of all 5 of 5 sequences from another bucket. Results are 
then sorted and displayed based on the computed hypergeometric-distribution 

10 statistic as shown in Figure 2 at step ST222. A number of standard 
algorithmic and bioinformatic optimizations can be made to improve system 
performance, such as but not limited to one or more of the following: pre- 
computing all the biomolecule relationships and using a look-up table to 
determine biomolecule identity or similarity, storing the subset of buckets with 

15 a significant number of matches in an associate array, and limiting the 
statistical computation to that subset. 

The problem of correlating a given query set with a target database is 
addressed. The methods and data structures disclosed herein can be readily 
implemented and employed in a range of applications. Additionally, the 

20 methods are able to tolerate small numbers of "contaminant" sequences in a 
bucket without significantly degraded performance. 
II.A. Construction of Target Database 

One property of the method is the generality of its application. Given 
any source of informative content about a particular set of biomolecules that 

25 share one or more common properties, the methods can create appropriate 
buckets and add them to an iterative, ever-expanding, and evolving target 
database. A target database thus comprises various classifications of 
biomolecules (e.g., genes and gene products) into collections, also known as 
"buckets", of entities having one or more common properties. 

30 A target database can be constructed. For example, as shown in 

Figure 4, a target database can be constructed by identifying a source of 
informative content (box 402 in Figure 4), arranging the informative content 
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into a set of named buckets (box 404 in Figure 4) wherein the members of 
each bucket share a common property, gathering the names of the buckets 
and a list of the biomolecules present in each bucket; and creating and 
loading into a database several data fields containing data representing the 
5 set of buckets, the list of biomolecules present in each bucket, a description 
for each biomolecule present in each bucket; an organism source for the 
biomolecule; and the user who inputted the information (see e.g., boxes 406 - 
412 in Figure 4). This data can be present as a data structure having multiple 
data fields and stored on a computer-readable medium, as is generally 

10 referred to as 400 in Figure 4. Interconnections between the data fields are 
schematically depicted by dashed lines in Figure 4. 

Referring now to boxes 408 and 410 in Figure 4, in one embodiment, a 
bucket can comprise a unique name describing its contents (e.g., "kinases"), a 
list of its members, and the nucleic acid and/or amino acid sequences for 

15 each of its members. A nucleic acid or amino acid sequence can be stored in 
one example as a file containing the nucleic acid or amino acid sequence 
itself. In another example, a nucleic acid or amino acid sequence can be 
stored as an identification number or accession number instead of the 
sequence itself, wherein the identification number or accession number allows 

20 the corresponding nucleic acid or amino acid sequence to be accessed as 
needed from a public or private database. For example, if the human 
erythropoietin gene or gene product is a member of a bucket, ( it could be 
stored in that bucket as the entire nucleotide or amino acid sequence of the 
human erythropoietin gene or protein, respectively. 

25 Continuing with boxes 408 and 410 in Figure 4, alternatively, the 

appropriate NCBI or SwissProt accession number can be stored instead. A 
variety of biological information including nucleotide and peptide sequence 
information is available from public databases provided, for example, by the 
National Center for Biotechnology Information (NCBI) located at the United 

30 States National Library of Medicine (NLM). The NCBI is located on the World 
. Wide Web at uniform resource locator (URL) "http://www.ncbi.nlm.nih.gov/", 
and the NLM is located on the World Wide Web at URL 
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"http^/www.nlm.nih.gov/". The NCBI website provides access to a number of 
scientific database resources including: GenBank, PubMed, Genomes, 
LocusLink, Online Mendelian Inheritance in Man (OMIM), Proteins, and 
Structures. A common interface to the polypeptide and polynucleotide 
5 databases is referred to as Entrez which can be accessed from the NCBI 
website on the World Wide Web at URL "http://www.ncbi.nlm.nih.gov/Entrez/" 
or through the LocusLink website. For the human erythropoietin gene and 
protein, for example, these accession numbers are AF202306 and P01588, 
respectively. In yet another example, the sequences can also be entered into, 

10 and subsequently retrieved from, a separate BLAST-formatted database. 
Each bucket entry can also contain a term describing the organism from 
which the reference sequence was derived (e.g., box 406 in Figure 4). Each 
bucket entry can also contain additional information, such as standard 
nomenclature for the gene or protein represented by the bucket entry. 

15 Continuing with Figure 4, in another embodiment, provided is the 

incorporation of a user-created query set into a target database. In this 
example, each set of nucleic acid or amino acid sequences that is submitted 
as a query can itself become a new bucket. In this example, the identity of 
each user, box 412 in Figure 4, can be tracked and the user required to enter 

20 the appropriate data into a common gateway interface (CGI) script-generated 
webpage. These "user buckets" can be treated in the same way as any other 
bucket. User buckets will likely be of varying quality, with some user buckets 
resulting from a thorough data analysis while other user buckets might be 
exploratory queries. 

25 The addition of user buckets can result in an enhancement in a given 

target database. For example, it is possible to add any (or all) gene clusters, 
dose-response or time-course gene sets, and lists of genes with altered 
expression derived from any experiment to a target database. Such additions 
can be made available to an entire project, group, site, or a corporate entity. 

30 Further, by identifying the user responsible for adding a specific user bucket, 
(e.g., by using bucket source identifiers as discussed hereinbelow), any user 
who finds that his or her query set is similar to that of another user will be able 
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to immediately recognize this event and notify the other user. Thus, 
communication of experimental results (e.g., results related to the implication 
of genes or gene products in different disease conditions) can be enhanced. 
Continuing with Figure 4, as the collection of buckets in a target 

5 database grows, it can be advantageous to define a "bucket source" 414 that 
describes the origin of each bucket. This can be desirable because two or 
more sources can often define buckets with exactly the same names, but with 
varying degrees of overlap in the sequence(s) that form the buckets. By 
including the source in a bucket name to form a "bucket source" identifier, the 

10 uniqueness of bucket names is assured. A further advantage of defining 
bucket sources is that it also facilitates defining broad categories of buckets, 
such as "pathway" or "function" buckets. This can be useful for helping to sort 
the output results or allowing users to choose and employ category types {i.e., 
buckets) that are most interesting or relevant to their work. 

15 An ad hoc rating system for relative ranking of the quality of each 

bucket source is optionally employed. In this rating system, reliable human- 
annotated sources (e.g., SwissProt accessible via the World Wide Web at 
http://us.expasy.org/sprot/) can receive a higher rank than a set generated by 
an automated computational procedure. • 

20 Continuing with box 414 in Figure 4, each bucket source can also have 

an associated file or URL comprising the raw data from which a given bucket 
was created, as well as a Perl script that parses the data and actually creates 
the bucket files. Since many data sources can be derived from public (e.g., 
SwissProt, TrEMBL, NCBI) or private databases (e.g., intra-corporation) that 

25 are continually changing, automated scripts can be employed for updating a 
target database collection periodically. There can also be some sets of 
buckets that are created once and need not be updated further. All buckets in 
the target database are stored in a relational database. A relational database 
enables rapid retrieval of data on any given biomolecule or bucket through the 

30 use of indexing. 

II. B. Comparison of a Query Set With a Target Database 
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The comparing of a query set (e.g., a user-defined set of nucleic acid 
or amino acid sequences) with a target database is disclosed, as is scoring 
and ranking the matches, and reporting the results. 



5 II.B.1. Searching a Target Database 

Pre-computed relationships of identity or similarity between 
biomolecules from other sources can be used. The identity relationships can 
be based on equivalence of accessions, identifiers, or names of genes and 
proteins from data sources such as NCBI's LocusLink, Swissprot, or HUGO. 

10 Thus, any member of the query set with a name, accession, identifier, or 
sequence identical to one in the target database can be considered a match. 
This identity relationship can be determined by the use of associative arrays, 
string matching, or regular expressions. More domain specific techniques 
might be applied for biomolecule sequences, such as BLAST (Altschul et a/., 

15 1990) or dynamic programming. A database of these pre-computed 
relationships or a method for computing these relationships that determines 
the identity or similarity of a member of the query set to that of a reference 
biomolecules can be employed. 

When all of the sequences comprising a query set and a target 

20 database comprise nucleic acid and/or protein sequences, the BLAST 
algorithms (Altschul et a/., 1990) can be employed to rapidly perform pairwise 
nucleic acid-nucleic acid, protein-protein, or nucleic acid-protein comparisons 
between each member of the query set and each member of a target 
database. In one embodiment, stringent BLAST parameters can be employed 

25 to enforce a strict matching criterion, thereby reducing the comparison to a 
binary response (i.e. match/no match) for each sequence pair. Stringent 
BLAST parameters can include, but are not limited to, parameters that require 
that in order for a match to be scored, two sequences must be sufficiently 
identical (e.g. 95%, 96%, 97%, 98% or greater) and the match region must be 

30 sufficiently long (e.g. 100 or more residues, or encompassing the entire length 
of any biomolecule of less than 100 residues in length). In this regard, it is 
noted that each target database match is not only a match to a specific 
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sequence, but also a match to the bucket(s) of which the sequence is a 
member. 

BLAST is one approach to identifying a degree of similarity between 
two or more sequences. Software for performing BLAST analyses is publicly 

5 available through the National Center for Biotechnology Information (NCBI: 
http://www.ncbi.nlm.nih.gov/), and also can be licensed from Washington 
University, St. Louis, Missouri, United States of America 
(http://blast.wustl.edu). 

The basic BLAST algorithm involves first identifying high scoring 

10 sequence pairs (HSPs) by identifying short words of length W in a query 
sequence, which either match or satisfy some positive-valued threshold score 
T when aligned with a word of the same length in a database sequence. T is 
referred to as the neighborhood word score threshold. These initial 
neighborhood word hits act as seeds for initiating searches to find longer 

15 HSPs containing them. The word hits are then extended in both directions 
along each sequence until the cumulative alignment score falls a 
predetermined value below the maximum achieved score. Cumulative scores 
are calculated using, for nucleic acid sequences, the parameters M (reward 
score for a pair of matching residues; always > 0) and N (penalty score for 

20 mismatching residues; always < 0). For amino acid sequences, a scoring 
matrix is used to calculate the cumulative score. Extension of the word hits in 
each direction are halted when the cumulative alignment score decreases by 
the quantity Xfrom its maximum achieved value, the cumulative score goes to 
zero or below due to the accumulation of one or more negative-scoring 

25 residue alignments, or the end of either sequence is reached. The BLAST 
algorithm parameters W, f, and X determine the sensitivity and speed of the 
alignment. For BLASTN searches, parameter settings such as 'W=13; 
hitdist=28; M=1; A/=-2; Q=1; R=1; X=6; gapw=20" from Washington University 
BLAST (WU-BLAST: http://blast.wustLedu/blast/TO-FLY.html#blastn) can be 

30 employed. For protein searches, parameter settings of "W=4; 7=1000; 
matrix=PAM10; £=1e-10" to determine identities can be used. 
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ln addition to calculating percent sequence identity, the BLAST 
algorithm also performs a statistical analysis of the similarity between two 
sequences. See, e.g., Karlin and Altschul, 1993. One measure of similarity 
provided by the BLAST algorithm is the smallest sum probability (P(N)), which 
5 provides an indication of the probability by which a match .between two nucleic 
acid or amino acid sequences would occur by chance. For example, a test 
nucleic acid sequence is considered similar to a reference sequence if the 
smallest sum probability in a comparison of the test nucleic acid sequence to 
the reference nucleic acid sequence is less than about 0.1, in another 
10 example less than about 0.01, and in still another example less than about 
0.001. 

Percent similarity of a DNA or peptide sequence can also be 
determined, for example, by comparing sequence information using the GAP 
computer program, available from the University of Wisconsin Genetics 

15 Computer Group (now part of Accelrys Inc., San Diego, California, United 
States of America). The GAP program utilizes the alignment method of 
Needleman and Wunsch (1970), as revised by Smith and Waterman (1981). 
Briefly, the GAP program defines similarity as the number of aligned symbols 
(i.e., nucleotides or amino acids) that are similar, divided by the total number 

20 of symbols in the shorter of the two sequences. See, e.g., Schwartz and 
Dayhoff, 1979, pp. 357-358, Gribskov and Burgess, 1986. 



II.B.2. Counting Matches 

In another aspect, guidelines are provided for counting matches. For 
25 example, when a candidate biomolecule of a query set matches a reference 
biomolecule of a target database (i.e. meets or exceeds a user-defined 
stringency requirement), a match is counted. When counting matches 
between sets, the following guidelines for counting matches can be employed: 

a) each query set candidate biomolecule can be considered to match 
30 at most one reference biomolecule in any given bucket; 

b) each query set candidate biomolecule can possess a match in one 
or more different buckets; and 



WO 2005/003308 



PCT/US2004/019932 



-29- 

c) once a candidate biomolecule in the query set matches a specific 
bucket reference biomolecule in the target database, any 
subsequent matches of that same candidate biomolecule to other 
reference biomolecules, or of other candidate biomolecules to the 
5 same reference biomolecule in that bucket, do not increase the 

match count for the bucket. 
The third guideline ensures that for a query set with Q members and a 
bucket with B members, the two cannot share more matches than the 
minimum of B and Q. A result of a counting procedure is a list of all the 
10 buckets in a target database that have one or more matches to a given query 
set. 



II.B.3. Statistical Significance of the Number of Matches 

In another aspect, the number of matches between a member of a 

15 query set and a bucket of a target database identified and counted as 
described herein can be analyzed to determine the statistical significance of 
the match. That is, the number of matches can be analyzed to determine, 
generally speaking, the likelihood that the number of matches is due to 
random coincidence, as opposed to a true property in common between the 

20 query set and the bucket of the target database. 

In general, the significance of a match will depend on the size of the 
query set, the size of each target database bucket that matched, the number 
of matches, and the total size of the relevant universe of all characterized 
sequences (approximated by the number of unique biomolecules in the 

25 reference collection). By way of example, the significance of a match can be 
modeled on the basis of a hypergeometric distribution as follows. If Q is the 
size of the query set, B is the size of a particular target database bucket that 
matched the query set, k is the number of matches between Q and B, and G 
is the size of the relevant universe, then the probability of exactly k matches is 

30 given by a hypergeometric distribution, which is defined as: 



P(Q D B = k) = C(B 9 k)C(G -B,Q-k)l C(G, Q) 
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where C(x,y) is the binomial coefficient: x!/[y!(x-y)!] and x! indicates factorial 
(the product of all integers from 1 to x). The P-value is indicated by the tail of 
the distribution: 

5 P(Q[\B :> k) = ^{c(B 9 k)C(G-B,Q- k)/C(G,Q)) 

The parameter G can be fixed as a constant for all computations. 

Although draft forms of the sequence of the human genome are 
available to the public for searching (See http://www.ncbi.nlm.nih.gov 
/genome/guide/human/), there is still no agreement on the number of genes or 

10 proteins it encodes. See e.g., Smaglik (2000), and compare Lander et al. 
(2001) estimating 30,000-40,000 protein-coding genes to Venter et al. (2001) 
estimating 38,000 genes. Although current estimates of the number of human 
genes range upwards from 20,000, many genes have not yet been 
characterized while others are likely incorrectly or incompletely characterized. 

15 There is no doubt that some genes have not yet been identified, so the true 
number of genes in the human genome is likely to be uncertain for many 
years to come. Therefore, any estimate made with respect to genome size is, 
to some extent, arbitrary. 

An aspect, therefore, pertains to the characterization of the number of 

20 genes comprising the human genome as a number reflecting how many 
human genes have been identified, annotated, or otherwise classified. 
Regardless, the specific value for the genome size has no impact upon the 
rank order of the buckets that are reported as significant matches. This 
degree of uncertainty in the size of the genome only affects the cutoff level for 

25 statistical significance. Thus, the relative ordering of the buckets is unaffected 
by any assumptions made concerning the size of the genome. It is also 
possible to compute an effective size for the genome of any organism by 
counting up all the unique sequences from that organism that have been 
partitioned into one or more buckets. Similarly, one could restrict the genome 

30 size to the number of probe sets (or number of unique genes) available on a 
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specific DNA microanray or chip, for purposes of analyzing experimental data 
from RNA expression studies. 

In one embodiment, the results of comparing a query set to a target 
database can be presented as a list of buckets ranked by p-value, and can be 
5 bounded by a predefined statistical cutoff. In one embodiment, for each of the 
buckets in the results list, a hyperlink can be incorporated in an output display 
that takes the user to a summary page. The summary page can be 
configured to show which query set sequences matched which bucket 
elements, as well as which bucket elements had no matches in the query set. 
10 One or more additional hyperlinks can also be included. These hyperlinks 
include, but are not limited to links to a database entry for each query set 
sequence (such as a link to the entry in SwissProt, NCBI, or a private 
database). 

15 II.B.4. Representative Steps 

The following section describes an embodiment of the method. The 
section generally describes a series of steps that can be performed when 
practicing the disclosed method. The following steps describe only one 
example. Variations on the disclosed method will be apparent to those of 

20 ordinary skill in the art, upon consideration of the present disclosure, and are 
encompassed by the appended claims. Reference is also made to Figure 2, 
where a method is referred to generally at 200. 

First, as shown at step ST202 and/or ST204a and ST204b in Figure 2, 
a query set comprising one or more candidate biomolecules is inputted to a 

25 computer that will run an analysis. A query set can comprise, for example, 
one or more sequences known or suspected to be located in the same 
genetic region. Alternatively, a query set can comprise an amino acid 
sequence of a protein known or suspected to be involved in a given biological 
pathway or complex, or can comprise a set of nucleotide or protein sequences 

30 which result from a biological experiment, such as gene or protein abundance 
changes, protein-protein interactions, etc. For convenience, sequences are 
inputted in the standard FASTA format. See Pearson, 1988 and Pearson, 
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1990. If the sequences of a query set are not in FASTA format, they can be 
converted to FASTA format. Additionally, the inputting can comprise entering 
accession identifiers and retrieving FASTA formatted sequences based on the 
identifiers, as depicted in steps ST204a and ST204b in Figure 2. 

5 Next, as shown in steps ST206 and ST208 in Figure 2, a sequence I of 

one or more candidate biomolecules of a query set is compared with a 
sequence of one or more reference biomolecules of a target database, the 
one or more reference biomolecules of the target database grouped into one 
or more buckets J. The comparison can be made using a matching of 

10 equivalent biomolecules names, sequences or accession, or by BLAST-based 
identity/similarity search based on sequence. Such a search can employ the 
algorithms of the BLAST method. Alternatively, the search can employ 
modified BLAST algorithms. The selection of the search algorithms to be 
employed can be made based upon consideration of the sequences and the 

15 target database composition. A target database can be generated as 
disclosed herein. Buckets can also be generated as disclosed herein. 

After a search has been performed, as shown in step ST208 of Figure 
2, a number of matches between each candidate biomolecule, e.g., sequence 
I, of the query set and each reference biomolecule in the target database are 

20 counted. This operation can generate a list detailing which candidate 
biomolecules of a query set matched which buckets, e.g., bucket J, (and 
which reference biomolecules) of a target database. Guidelines for counting 
matches are provided herein. The iterative capability of the present method is 
shown in steps ST210, ST212, ST214, ST216, and ST218, wherein the 

25 search continues to an additional bucket J+1 . 

Following a search, as shown in step ST220 of Figure 2, one or more 
statistics (such as hypergeometric statistics or hypergeometric statistics 
including empirical correction multiple hypothesis testing) for each bucket 
match can be computed. Such a computation can account for the genome 

30 size G, and the query set size Q, and can be based on bucket size B and 
number of hits k. By performing a very stringent BLAST search, it is possible 
to assert whether or not a sequence is present in any given bucket, and 
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should therefore be considered as possessing that biological property. By 
performing a statistical analysis, it can be determined how likely it would be 
for a similar result to have been obtained by chance, if choosing at random 
the same number of sequences in the input set This enables the results to 

5 be ranked, with those properties shared by all (or a large subset) of the query 
set to receive greater priority than those properties that only occur in an 
individual sequence from the query set. 

Standard cut-offs for p-values (such as 0.05, 0.01, 0.001) can be used 
to guide significance. These p-values can be corrected for multiple 

10 hypothesis testing using a suitable approach, such as but not limited to one or 
more of a conservative Bonferroni correction (which multiplies these values by 
the number of hypotheses tested equal to the number of buckets for this 
embodiment) and computing an empirical p-value based in simulations with 
random input sets. This empirical p-value can be obtained by using multiple 

15 random input sets of genes and computing the number of times any bucket is 
observed below a certain statistic. For example, the algorithm can be 
simulated 1000 times on random input sets of genes (each set with 50 
members). The distribution of the best observed hypergeometric statistic from 
each of those 1000 computations can be plotted, and a statistic chosen, such 

20 that only 50 of the 1000 simulations have a statistic as good. This effectively 
gives the statistic that represents an empirical p-value of 0.05 for query sets of 
size 50. This can be repeated for query sets of varying sizes. 

The results of the statistical operation can then optionally be sorted by 
increasing or decreasing significance, as shown in step ST222 of Figure 2. 

25 The results of the operation can then be displayed to a user. Convenient 
display formats can include an output webpage. When results are displayed 
on a webpage, the results can be accompanied by hyperlinks to further details 
of the search, to the match, to the target database, and/or to the query set 
members. 

30 



WO 2005/003308 



PCT/US2004/019932 



-34« 

ML Genomic Region Analysis 

The methods disclosed herein can be employed to identify a property 
common to a set of candidate biomolecules from one genomic region that 
form a query set and a set of reference biomolecules that form one or more 

5 buckets of a target database. However, the present method is not limited to a 
comparison of a query set comprising a single set of candidate biomolecules 
and a target database. As described in the following sections, one 
embodiment of the method can be employed to identify a property common to 
a query set comprising two or more region sets and a target database. 

1 0 Representative steps are as follows: 

a. providing a query set describing two or more region sets, each 
region set comprising one or more candidate sequences extracted 
from a genomic region; 

b. comparing the query set with target database sequences describing 
15 one or more reference biomolecule sequences, the target database 

sequences grouped into one or more buckets, and wherein the one 
or more reference biomolecules of each bucket share a common 
property; 

c. counting a number of matches between each query set and each 
20 bucket of the target database; and 

d. statistically analyzing each match, wherein the presence of a 
statistically significant match identifies a relationship between the 
query set and the bucket of the target database. 

A non-limiting example of this embodiment can be described in the 
25 context of a disease gene association analysis, and is referred to generally at 
300 in Figure 3. As shown in step ST302 in Figure 3, a given set of genomic 
regions known or suspected to be associated with a particular disease or 
general disease category is first identified. As shown in steps ST304 and 
ST306, for each such region, endpoints are determined and a set containing 
30 at least some and preferably all of the known and predicted genes that lie 
within the region is created. This set is known as a tt region set". Multiple 
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region sets can be accommodated. Region sets can be combined to form a 
query set. 

A query set, which comprises two or more region sets, is then 
compared, region set by region set, with a target database, at step ST308 in 
5 Figure 3. The comparison can optionally be made by one of employing either 
the equivalency of name, identifier or accession, and by using BLAST or 
BLAST-based algorithm(s) on the sequences of the biomolecules. For each 
region set in a query set, if one of the candidate sequences of a region set 
matches a reference sequence in the target database, the matching 

1 0 sequence(s) is scored as "present" for that region. 

Continuing with step ST308 in Figure 3, each candidate sequence of a 
query set is sequentially compared with the reference sequences of the target 
database to generate a list of candidate sequences in the query set (which 
can be sorted by region set) that are present in the target database. The 

1 5 query set sequences that match a reference sequence in a target database 
can be sorted by target database set(s) {i.e., buckets) that contain at least one 
match from a specified number of different region sets. This process 
generates a list of buckets found in one or more region sets. 

As shown at steps ST310, ST312, and ST314 in Figure 3, the 

20 statistical significance of a match between a query set sequence and a target 
database set can then be calculated. The method can also be adapted to 
allow the results to be sorted and displayed on the basis of one or more 
criteria. The incorporated statistical analysis offers a step for ensuring that 
any observed result (e.g., a match between a query set sequence and a 

25 target database sequence) cannot be explained solely by random chance. In 
one example of such a statistical analysis, simulations are employed to 
randomly choose a set of genomic regions with similar gene numbers to the 
input data, to compare the simulation data to the sequences of a target 
database, and to score any matches observed. Subsequently, another 

30 random set of genomic regions is chosen and the process is repeated until a 
predetermined number of iterations or replicates has been performed. In one 
example, and as depicted in step ST316 of Figure 3, iterations are done for 
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1000 random region sets. Then, the results of the simulation are compared 
with data obtained from an actual query set, as shown in step ST318 in Figure 
3. Actual data set matches that rank statistically highly in the simulation can 
be considered to be potential false positives and can be discarded as not 
5 indicative of a meaningful match. The results of the statistical analysis can 
then be displayed, as shown in step ST320 in Figure 3. For example, for a 
bucket appearing with statistic of 0.01 in the actual analysis, if in the 1000 
simulations that bucket appears with that statistic or better only 50 times, then 
this bucket is assigned an empirical p-value of 0.05. This would be less 
10 significant than the p-value based on the theoretical statistic alone. 

IMA Searching a Target Database 

The steps summarized immediately above will now be discussed in 
detail. The identification of one or more properties of a candidate sequence of 

15 a query set can be achieved by searching a target database of reference 
sequences that have been grouped into buckets representing groups of 
sequences that have the same properties. Such a search can follow another 
experiment, the results of which can form elements of a query set. Some 
technologies and experiments generate a powerset of genes, for example s 

20 =(SiS 2 . . .S n ). A subsequent goal is then to find a property P such that there is 
at least one gene in a significant number of sets S k that has property P. The 
sets Si...S n have no pairwise intersection (e.g., non-overlapping genomic 
regions). Thus if biological pathways are considered as potential properties, 
then the goal might be to find a pathway that threads or connects these sets 

25 of genes. 

Consider a linkage analysis experiment, where genetic markers have 
been genotyped in both a disease and a normal population and log of the 
odds ratio (LOD) scores have been obtained across the human genome. For 
most common diseases, multiple linkage peaks are observed. The presence 
30 of multiple linkage peaks can be explained as: 

• none of the linkage peaks are real; OR 
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• some of the higher linkage peaks are real, and each peak has one 
or more genes in the corresponding genomic region, mutations in 
which explain the observed linkage 



5 If the latter proposition is true, then one of the following cases must 

apply. 

• these mutated genes are independent (they do not share any 
property, such as a pathway) 

OR 

10 • (H) a subset of these genes is related by a property, such as a 

pathway, and mutations in a subsection of this pathway all 
contribute to the same phenotype/disease 

MLB. Counting Matches 

15 Exploration of known properties that might explain hypothesis H above 

can be accomplished using a genomic region analysis embodiment. Consider 
n genomic regions with corresponding region sets of genes: Si,S 2 ...S n . If an 
additional set of genes R is added to (S 1f S 2 ...S n ), where R contains all 
remaining genes not in Si...S n , then the superset (Si...S n , R) can be 

20 considered a partition on the human genome. Thus, consider the data 
superset s = (Si...S n , R) and let the number of genes with property P 
overlapping with these sets be (Pji...P Jn , Pj r ). The probability of this event (or 
partition j) is given by the multivariate form of the hypergeometric distribution 
(sampling without replacement): 

25 Prevent j) = C(R,P Jr )f[C(S k ,P Jk )/C(Genome,\ P |) 

where C() is the binomial coefficient, and |S| is the cardinality {i.e., the number 
of members) of the set S. The probability of seeing this by chance can be 
estimated by summing the above term over all events -j that would be 
considered significant, for examples, events that have at 3 or more P jk greater 
30 than 0. 
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ln certain applications, it might only be important that the region set S 
have at least one biomolecule in common (or identical) to that with property P. 
It might not add any more evidence if it has two or more molecules with 
property P. In such cases, computing an exact significance (p-value) 
5 becomes a difficult task, and Monte-Carlo techniques can be used to acquire 
estimates as discussed in the next section 

III.C. Statistical Significance of Matches 

10 Statistical measures make assumptions of independence among set 

members that do not completely hold for biological sequences. Thus, the 
significance of any result is assessed using negative controls. True negative 
controls are hard to obtain, as that would require knowing that a certain 
powerset of biological sequences shares no property in a significant measure. 

15 A solution is to generate multiple random sets and use simulations to compute 
the background frequency of a property. A similar approach is adopted here. 
A plurality of replicates of a query set is constructed. These replicates are 
matched to the query set in the sense that the number of genes in each 
replicate is equal to the number of genes in each set of the powerset and as 

20 far as possible arises from a similar bioinformatics process. The replicates 
are then modeled at random chromosomal locations to form a random 
location data set. The random location data set is then processed using the 
same method steps described above. For example, if the original powerset 
represented linkage regions, then each random set would be a set of 

25 contiguously ordered genes from a single chromosome, and a random set of 
genes from a contiguous region of the genome can be generated. For each 
property, the number of times that property is observed in the random 
powerset is counted with P(event) equal to or lower than that observed in the 
actual powerset. This provides a simulation-based or empirical p-value. 

30 A similar approach can also be used in the analysis of time-series data, 

such as data gathered from microarray expression experiments over time. 
For example, some experiments produce a list of genes with significantly 
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perturbed expression at the earliest time point, and various other sets that 
experienced expression changes at successively later time points. Some 
pathway or process that connects this set of measurements can also be 
identified. As an example, assume that there were three time-points 
5 measured - (E)arly, (l)ntermediate, and (L)ate. For each time point there 
would be an associated set of genes (E|, l it Li) whose expression levels had 
changed relative to the control (time=0). If these sets are non-overlapping, 
one can apply the method to discover any processes that contain one or more 
genes that are present in each timepoint set, thus forming a hypothesis as to 

10 the pathway and the causal steps involved in the experimental process. 
Statistical corrections are employed to handle the case where the sets are not 
completely disjoint. 

By way of additional example, schizophrenia is a multifactorial disease. 
A number of linkage studies have been published implicating the following 

15 chromosomal regions: 1q21-22, 1q32-42, 6p24-22, 8p21, 10p14, 13q32, 
18p11, and 22q11-13. Blouin et a/., 1998; Berrettini, 2000; Straub et a/., 
1995; Brzustowicz et a/., 2000; Ekelund et a/., 2001 . For the chromosome 1q 
region, conflicting evidence also exists. Levinson et a/., 2002. Given suitable 
markers or other methods to determine the physical boundaries of each 

20 region, one can extract the set of known and predicted genes within each 
such region. The genome region analysis' embodiment is then used to probe 
for pathways or other biological processes that have components in some or 
all of the linkage regions. Simulations can also be performed to repeatedly 
generate randomly located chromosomal regions of comparable size and 

25 gene content to assess whether the results occur frequently by chance alone. 
The findings are then used as hypotheses for guiding experimental studies. 

Examples 

The following Examples have been included to illustrate certain 
30 embodiments. Certain aspects of the following Examples are described in 
terms of techniques and procedures found or contemplated to work well in the 
practice of the embodiments. These Examples are exemplified through the 
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use of standard practices of applicants. In light of the present disclosure and 
the general level of skill in the art, those of skill will appreciate that the 
following Examples are intended to be exemplary only and that numerous 
changes, modifications, and alterations can be employed without departing 
5 from the spirit and scope of the present disclosure. 
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Example 1 
Pseudocode 

• If input not FASTA format, read accessions and get FASTA sequences, 

• Compare input sequences against entire bucket database (use BLAST- 
5 based identity search or simply accession ID lookup), 

• For each input sequence, count number of matches to each bucket in 
database, 

• Given the genome-size G, and the query set-size Q, compute 
hypergeometric statistic for each bucket possessing matches, based 

10 on bucket-size B and number of hits k. 

• Sort the results list by decreasing significance and output webpage 
with results and hyperlinks to further details. 



Example 2 

15 Analysis of Genes Regulated by E2F 1 

Stanelle reported 29 genes as being regulated by the transcription 
factor E2Fi. Stanelle et a/., 2002. The authors divided this set of genes into 
five categories: cell cycle, apoptosis, cancer-related, E2Fi targets, and 
unknown. Submitting the same unordered list in an embodiment of the 

20 present method results in a ranked list of approximately 100 buckets 
significant at p £ 0.05. Presently there are approximately 80,000 buckets in 
the target database. These buckets have been created from a combination of 
publicly available databases and internal experimental results. These buckets 
cover many types of biological data including, but not limited to genomic 

25 location, diseases, tissue expression, functions, pathways, transcriptional 
regulation, families, domains, and literature abstracts. The most significant 
hits of this input set to the target database are shown in Table 1 . Some of the 
sources which appear are keywords and families from Swissprot, protein 
domains from EnsEMBL Interpro, human disease sets from OMIM, and sets 

30 derived from the Gene Ontology Consortium website and NCBI's LocusLink. 
This list includes several overlapping buckets related to each of the known 
categories supplied by the authors, with cyclin C (cell-cycle) determined to be 
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the most significant bucket In addition to confirming the authors' 
classifications, more specific links, such as associations with MAPKKK 
signaling and multiple myeloma, were also uncovered. See Table 1 . 



5 Table 1 

Results of Classifying the Set of Genes Regulated bv 
Transcription Factor E2F i 





Bucket Name 


Bucket Source 


Bucket Size 


Obs 


Exfi 


D value 


10 


CyclinC 


EnsEMBL InterPro Domains 


12 


4 


0.046 


8.7e-08 




Apoptosis: BP 


LocusLInk GO 


227 


9 


0.88 


1.16-07 




Cell Death: BP 


LocusUnk GO 


236 


9 


0.91 


1.5e-07 




Death: BP 


LocusUnk GO 


236 


9 


0.91 


1.5e-07 




Induction Of Apoptosis: BP 


LocusUnk GO 


90 


6 


0.35 


9.6e-07 


15 


G1/S-Specific Cyclin: MF 


GOA 


9 


3 


0.035 


4.3e-06 




Cydln 


EnsEMBL InterPro Domains 


33 


4 


0.13 


6.86-06 




Cyclin 


SwissProt Keyword 


18 


3 


0.07 


4.1e-05 




E2F1 


OMIM 


18 


3 


0.07 


4.1e-05 




Cyclin 


SwissProt Family 


19 


3 


0.073 


4.8e-05 


20 


Cyclin: MF 


GOA 


19 


3 


0.073 


4.8e-05 




Apoptotic Mitochondria) 


LocusLink GO 


4 


2 


0.015 


8.6e-05 




Changes: BP 














Apoptotic Program: BP 


LocusLInk GO 


25 


3 


0.09 


0.0001 




Myeloma, Multiple 


OMIM 


5 


2 


0.019 


0.00014 


25 


Induction Of Apoptosis By 
Extracellular Signals: BP 


LocusLink GO 


30 


3 


0.12 


0.00020 




MAPKKK Cascade: BP 


LocusLink GO 


34 


3 


0.13 


0.00029 



30 Example 3 

Pseudocode Genomic Region Analysis Embodiment 
For each genomic region of interest, extract at least some and 
preferably all of the known genes contained therein. 
For each region set, compare each candidate sequence to the bucket 
collection (use BLAST-based identity search or simply accession ID 
lookup). 

For each bucket in the database, count number of region sets that 
contain at least one biomolecule in common with the bucket. 



35 
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• Choose some constant M <> number of regions, and report all buckets 
that had hits to at least M regions. 

• Use multivariate form of hypergeometric distribution to assess 
significance of these buckets. 

5 • Given the number of regions and number of genes in each region, 
construct 1000 replicates of the region set (same number of regions 
and same number of genes per region), but placing the simulated 
regions at random chromosomal locations. 

• Process this random data set in the same way as the real data, and 
10 note how many times (out of 1000 replicates) that each bucket was 

found with at least as strong statistical support as it had in the original 
data set. 

• Display the results, ranked by decreasing statistical support in the 
original data. 

15 
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20 It will be understood that various details can be changed without 

departing from the scope of the disclosure. Furthermore, the foregoing 
description is for the purpose of illustration only, and not for the purpose of 
limitation. Indeed, variations on the present disclosure will be apparent to 
those of ordinary skill in the art, upon consideration of the present disclosure, 

25 and are encompassed by the appended claims. 



